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ABSTRACT 

We evaluate the coorbital corotation torque on a planet on a fixed circular 
orbit embedded in a viscous protoplanetary disk, for the case of a steady flow in 
the planet frame. This torque can be evaluated just from the flow properties at 
the separatrix between the librating (horseshoe) and circulating streamlines. A 
stationary solution is searched for the flow in the librating region. When used 
to evaluate the torque exerted by the circulating material of the outer and inner 
disk on the trapped material of the librating region, this solution leads to an 
expression of the coorbital corotation torque in agreement with previous esti- 
mates. An analytical expression is given for the corotation torque as a function 
of viscosity. Lastly, we show that additional terms in the torque expression can 
play a crucial role. In particular, they introduce a coupling with the disk den- 
sity proflle perturbation (the 'dip' which surrounds the planet) and add to the 
corotation torque a small, positive fraction of the one-sided Lindblad torque. As 
a consequence, the migration could well be directed outwards in very thin disks 
(aspect ratio smaller than a few percent). This 2D analysis is especially relevant 
for mildly embedded protoplanets (sub-Saturn sized objects). 

1. Introduction 

The discovery of extra-solar planets in the last few years renewed interest in the 
description of the tidal interaction of a planet with a gaseous Keplerian disk. Indeed, many 
planets, the so-called 'hot Jupiters', have been found orbiting very close to their host star, 
and their formation is unlikely to have occurred at such short distances from the star. 
It is thought that these planets have formed further out in the disk and have undergone 
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a significant orbital decay, or migration, as the result of the tidal interaction with the 
protoplanetary disk in which they were formed. Goldreich & Tremaine (1979, 1980) and 
Ward (1986, 1997) have shown that most of the tidal interaction between a protoplanet 
embedded in a protoplanetary disk can be accounted for by the so-called Lindblad torques, 
which correspond to the exchange of angular momentum between the planet and the disk at 
the position of the Lindblad Resonances (LR) of the planet, where the gravitationnal field 
of this latter excites in the disk spiral density waves which propagate away from the planet's 
orbit. The Lindblad torque exerted by the outer disk (at the Outer Lindblad Resonances, 
OLR) on the planet is negative, whereas the Lindblad torque exerted by the inner disk (at 
the Inner Lindblad Resonances, ILR) on the planet is positive. Ward (1997) has shown 
that the absolute value of the Outer Lindblad torque is intrinsicly higher than the Inner 
Lindblad torque; hence the net torque (which we shall refer to hereafter as the differential 
Lindblad torque) is negative and leads to an orbital decay of the planet. Migration can 
roughly speaking be split into two regimes: 

• The type I migration regime, in which the planet mass is small, and the disk response 
is linear, that is to say that the relative perturbed density in the planet's wake is 
a small fraction of unity. In this regime, the migration velocity is proportional to 
the planet mass, and is independent of the viscosity, as long as one only considers 
the Lindblad torque and neglects the corotation torque (Meyer- Vernet & Sicardy 
1987, Papaloizou & Lin 1984). The LR pile- up at high azimuthal wavenumber at 
±|/i (where h is the disk thickness) from the planet's orbit, and most of the torque 
acting on the planet comes from the pile-up region. The relative mismatch between 
the Outer and Inner Lindblad torques, which gives rise to migration, scales as the 
distance between the pile-up of the ILR and the pile-up of the OLR, and therefore 
scales as the disk aspect ratio h' = h/r. Since the one-sided Lindblad torque scales 
as h'~^ (see e.g. Ward 1997), the migration velocity scales as h'~'^. 
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• When the protoplanet mass is above a certain threshold (namely when its Hill radius 
is comparable to the disk thickness), the wake excited by the planet leads to a shock 
in the immediate vicinity of the excitation region, and therefore the wake is locally 
damped and gives its angular momentum to the disk. Hence the planet can open a 
gap, if the viscosity is below some threshold value (see e.g., Papaloizou & Lin 1984). 
In that case, the disk response is markedly non-linear, and most of the protoplanet 
Lindblad resonances fall in the gap and therefore cannot contribute to the planet-disk 
angular momentum exchange. This is the so-called type H migration regime. The 
migration rate slows down dramatically compared to the type I migration (Ward 
1997). Furthermore, the tidal truncation splits the disk in two parts and, if the 
surrounding disk mass is comparable to the planet mass, the planet is locked in the 
disk viscous evolution. In this case the migration time is of the order of the disk 
viscous timescale (Nelson et al. 2000). 

One of the problems raised by type I migration is that it is extremely fast. In 
particular, the build up of a single giant planet solid core, which has to reach a critical mass 
of about 15 (where is the earth mass) before rapid gas accretion begins (Pollack et 
al. 1996), is a slow process, and can be longer by more than one order of magnitude than 
its migration timescale all the way to the central star (Papaloizou & Terquem, 1999). The 
planet can also exchange angular momentum with its coorbiting material. This corresponds 
to the so-called corotation torque, which is proportional to the vortensity gradient, i.e., 
to the quantity d{B/T,)/dr, where S is the disk surface density and B the second Oort's 
constant (Goldreich & Tremaine, 1979). On the contrary of what happens at the Lindblad 
resonances, the flux of angular momentum at corotation is not carried away by waves but 
accumulates there. Ward (1991, 1992) has shown how the corotation torque is due to the 
exchange of angular momentum between the planet and fluid elements as they jump on 
the other side of the planet orbit when they reach a hair-pin curve end of their horseshoe 
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streamline, and that the dependence of this torque on the vortensity gradient arises from 
an unequal mapping of surface elements between the outer and inner parts of a horseshoe 
orbit. Estimates of this torque, both analytical (Ward 1992) and numerical (Korycanski 
& Pollack 1993) show that it is at most comparable to the differential Lindblad torque. 
Furthermore, in an inviscid disk, the libration of the coorbital material tends to remove 
the vortensity gradient after a time which is of the order of the turn-over time of the 
outermost horseshoe orbit (which contribute most to the coorbital corotation torque), thus 
switching off the corotation torque after a transient episode. In a viscous disk, the situation 
is different. On one hand, as suggested by Ward (1991), if the viscous diffusion timescale 
over the horseshoe region width is smaller than the turnover time of the horseshoe orbits, 
the gradient of vortensity will not be removed. On the other hand, the material of the outer 
disk, when it reaches the coorbital region during its drift towards the center, exchanges 
once angular momentum with the planet as it gets sent to the inner disk. 

The purpose of this paper is to analyze such a situation. After introducing the problem 
set-up and the main notations in section 2, we describe the streamlines in stationary regime 
(section 3). We show that the introduction of a finite viscosity modifies the topology of the 
flow around the planet, but that there still exists a set of hbrating horseshoe streamhnes, 
which correspond to material trapped in the coorbital region, and which therefore does not 
participate in the disk accretion onto the central star. We then show in section 4, assuming 
that the flow is stationary in the frame corotating with the planet, how the coorbital 
corotation torque can be expressed using only the flow properties at the boundary of the 
librating region (which we shall call hereafter the separatrix, as it separates the librating 
from the circulating streamlines). This simplifies greatly the torque expression and 'hides' 
any additional complexity arising from the flow properties interior to the separatrix, for 
instance, the "tadpole" streamlines around the Lagrange points L4 and L5. 
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In Appendix A, we give a simple steady state solution for the perturbed density in 
the coorbital region which neglects pressure effects. When used to evaluate the torque 
expression given in section 4, we recover the expression given by Ward (1991). We then 
give an analytical expression of the corotation torque as a function of viscosity, which 
exhibits two different regimes, depending on the ratio of the turnover to viscous timescales 
of the horseshoe region. A third regime is also mentioned, which corresponds to a cut-off 
at high viscosity when some of the streamlines originating from the outer disk arrive at the 
inner disk without having passed by the planet, i.e., when the viscous drift time over the 
horseshoe region is shorter than the horseshoe turnover time. 

In section 6, we show that the full corotation torque expression as given in section 4 
exhibits additional terms, which have potentially important effects. Indeed, they correspond 
to a small positive fraction of the one-sided Lindblad torque, and this fraction turns out to 
be independent both of viscosity and disk thickness. Since the relative mismatch between 
Inner and Outer Lindblad torques scales as /i', this additional term can be larger in absolute 
value than the differential Lindblad torque in sufficiently thin disks. This term, which can 
be estimated only in order of magnitude, may therefore significantly slow down, or possibly 
reverse, type I migration in very thin disks. 



2. Problem set-up and notations 

We consider a non-accreting protoplanet on a fixed circular orbit of radius in a 
disk with uniform and constant aspect ratio h' = h/r and viscosity u. Since the planet 
is considered as non- accreting, its mass Mp has to be below a threshold of about 15 
(Papaloizou & Terquem, 1999). The disk surface density is denoted E, and the unperturbed 
disk surface density is uniform and is denoted Sqo- It is also the value far from the planet. 
The viscosity is also given by: u = aflh"^, where Q is the disk orbital frequency, and a 
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is a dimensionless parameter (Shakura & Sunyaev 1973). The planet orbital frequency is 
denoted Qp. A fluid element is labeled in the disk by its coordinates r and ^ in a frame 
corotating with the planet, where r is the distance to the primary and 9 the azimuth (the 
planet is at azimuth 0), counted rotation- wise. The perturbed azimuthal velocity in the 
presence of the planet is denoted v, and j — r^Q + rv is the speciflc angular momentum 
of a fluid element. Its value in the unperturbed disk is Jq = r^Q. We shall also make an 
extensive use of the distance to the orbit x — r — rp, and its non-dimensional counterpart 

3C — X j '^p • 

3. Streamline topology 
3.1. Inviscid case 

In an inviscid disk, the streamlines in the coorbital region of a planet on a flxed circular 
orbit have the aspect represented in Fig. la, lb and Ic. The shaded region represents the 
horseshoe streamlines. The streamlines in the outer and inner disk, which are represented 
as dotted hues in Fig. lb and Ic, are circulating, and closed in the case of a steady flow (as 
can be seen in Fig. la, where these streamlines appear as almost circular streamlines), since 
there is no radial drift (Lubow, 1990). The separatrix between the librating (horseshoe) 
and circulating (inner and outer disk) streamlines is represented as a thick solid line. The 
planet lies at r/r^ = 1 and ^ = (mod. 2n). The circulating streamlines always stay on 
the same side of the planet's orbit, as is better seen in Fig. Ic. One can see that these 
circulating streamlines get deflected as they pass by the planet. This corresponds to the 
excitation of epicyclic motion and the exchange of angular momentum with the disk, which 
is carried away by pressure effects and gives rise to the wake, which can be clearly seen in 
Fig. lb and Ic. The angular momentum transfer rate into this wake corresponds to the 
Lindblad torques. 
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Fig. 1. — The figures a to c represent the same situation: a 15 protoplanet embedded 
in a 4 % aspect ratio disk, once a steady state has been reached. The details of the run 
are unimportant here, as these figures are for illustration purpose only, and closely resemble 
other runs presented elsewhere (see e.g. Massct & Sncllgrove 2001, or Nelson et al. 2000). 
Figure d shows streamlines which fulfill the prescriptions 1 to 3 of section 3.2. The figures e 
and / represent the same situation, except for a 7r-shift in azimuth. The run they are 
extracted from differs from the run of figures a to c only by the viscosity, which has been set 
this time to a very high value. In figures d to /, the solid line in the outer and inner disk 
(outside of the shaded -librating- region) represents a unique streamline, which is also, as 
we are dealing with a steady flow, the path of a fluid element on its way to the central star. 
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3.2. Viscous case 

In order to have a first idea of the streamline topology we approximate the action of 
the planet on the flow with the following prescription: 

1. The planet acts on a fluid element at 6' = (mod. 2tt) only (i.e., at a conjunction). 

2. If the distance |a;| of the fluid element to the orbit when it is in conjunction with the 
planet, is smaller than some threshold value Xc (which corresponds to the half- width of 
the horseshoe region), then the fluid element is "reflected" with respect to the planet 
orbit and sent to {6, —x). Otherwise no action is taken on the fluid element. This is 
meant to mimic the U-turn of the fluid element at the end of the horsestreamlines in 
the horseshoe region. 

3. The fluid element velocity is assumed to be everywhere the velocity in the unperturbed 
disk. 



Obviously this oversimple prescription reproduces in an inviscid case the main 
properties outlined at the previous section, i.e., it gives rectangular librating streamlines 
in the horseshoe region (the width of which is 2xc), and circulating circular streamlines 
in the outer and inner disk. However there is no wake. In a viscous disk with uniform 
viscosity and surface density, the fluid element velocity in the corotating frame is given by 
(if we assume the radial extend of the coorbital region to be small compared to the planet 
distance) : 

3 v , • 3fl. 



which integrates as: 



X — — ^— and 9 — —-—x (1) 



0^1%' + 00 (2) 
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The streamlines are therefore arcs of parabola, and the sketch of the streamhnes which 
respect the prescription given above is represented in Fig. Id. If we follow the path of a 
fluid element along the solid streamline originating in the outer disk, we see that first this 
fluid element has a circulating trajectory and approaches progressively the planet coorbital 
zone. The distance between two successive intersections of the streamline with say the 
^ = axis increases steadily, since the angular frequency mismatch between the planet and 
the fluid element decreases. When the fluid element reaches the point A', it is not sent on 
the other side of the orbit, since its distance to the planet is still larger than Xc- However, 
along the arc A' A, it crosses the horseshoe region outer boundary, and when it reaches A, 
it is sent to B and therefore gives to the planet, during this close encounter, the positive 
amount of angular momentum that it loses. It then follows its path along BB', which by 
construction is symmetric to AA', therefore when it reaches B', it is out of the horseshoe 
region, and therefore it keeps circulating in the inner disk, and eventually gets accreted 
onto the primary. 

One can wonder how this schematic picture changes in a more realistic case. First it 
can be useful to notice that in a 2D steady flow defined on a compact domain, a streamline 
which does not intersect the domain boundaries is closed^. Therefore one has to expect a 
region of closed streamlines behind the planet, which enclose the fixed point on the orbit 
where the negative viscous torque cancels out the positive torque from the planet, which 
goes from to an arbitrarily large value as 9 goes from tt to 27r, and, if the viscosity is small 

^One reason for that is that the field T,u, where u is the velocity field, is divergence free, 
and can therefore be expressed as the curl of a field ipel, where ■0 is a scalar field and a 
vector orthogonal to the disk plane. Hence the function ip is conserved along the field lines 
of EtT, which are also the streamlines. The streamlines therefore appear as isocontours of 
the function ijj. 
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enough, this closed region should extend from one end to another of the horseshoe region. 
We represent in Fig. le and If a situation taken from a run similar to the run of Fig. la 
to Ic, except that the viscosity has been set to a very high value (corresponding to a ~ 0.1, 
which is probably too high for a protoplanetary disk, but has the merit that it clearly shows 
the streamhnes topology). We see on this figure that we have a closed streamhnes region 
which resembles the one obtained with the oversimple prescription, and that material flows 
from the outer disk to the inner disk. 

Therefore, in the case of a finite viscosity, we still have a set of closed librating 
streamlines in the coorbital region, enclosed by a closed separatrix, which corresponds to 
trapped material which will not be accreted onto the primary, but now the outer and inner 
disk communicate, and as material flows from the outer disk to the inner disk it exerts on 
the planet a positive torque. The smaller the viscosity, the more the librating region has a 
rectangular shape and resembles the horseshoe region of the inviscid case. 



4. Coorbital corotation torque expression 

The material enclosed within the separatrix is trapped and therefore its angular 
momentum is constant with time, since the planet is on a flxed circular orbit. Hence the 
net torque on the librating region L vanishes. This torque can be decomposed as the 
gravitational torque exerted by the planet Tp_^L and the viscous torque F^'* exerted by the 
inner and outer disk on L at the separatrix hence: 

Vp^L + Tf^Q (3) 

In what follows we make the following assumptions: 

1. the viscosity is small enough so that the separatrix distance to the orbit is almost 
constant when one stands far from the horseshoe U-turn regions; 
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2. the outer and inner separatrices lie at the same distance Xg from the orbit. This 
neglects in particular the mismatch between the orbit and the corotation, due to the 
fact that the disk is slightly sub-Keplerian (due to the partial support of the central 
object gravity by a radial pressure gradient). For an aspect ratio of a few percents, 
this mismatch is typically 10~^rp, while for the protoplanet masses we consider, ~ a 
few 10~^rp, and therefore assumption 2 is reasonable. However, when considering 
disks with aspect ratios > 0.1, and/or smaller mass planets, this mismatch can no 
longer be neglected. 

3. the perturbed quantities in the librating region are independent of 9, and thus are 
only function of x. 

Experience from numerical simulations shows that this last assumption is reasonable; for 
instance, the shallow 'dip' profile that is observed when the mass is sufficiently large and 
the viscosity is sufficiently low can be considered as independent of 9 with a very good 
approximation as long as one excludes the U-turn ends of the horseshoe region. We neglect 
the torque exerted by the inner and outer disk on L at the U-turn ends of the librating 
region, and we neglect the torque of the pressure forces on L. The coorbital corotation 
torque on the planet, which is the torque exerted by the disk material during the close 
encounters with the planet^, can therefore be expressed, using Eq. (3) as: 

Tc = ri^p + r« = rr + r^ (4) 

point of terminology is needed here: in an inviscid situation, the coorbital corotation 
torque simply corresponds to the exchange of angular momentum with the fluid elements 
which corotate, in average, with the planet, i.e. material strictly located at the corotation 
in the linear limit, and, in the case of a finite mass, the whole material of the horseshoe 
region (indeed, as a fluid element located in this region can never be in conjunction with 
the planet, when this latter has described N orbits in an inertial frame, the fluid element 
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where T^^p — —Tp^i is the gravitational torque exerted by the Ubrating region on the 
planet, Fr is the torque exerted on the planet by the circulating material when it flows 
from the outer disk to the inner disk, T'"g+ is the torque exerted by the outer disk on the 
librating region at the outer separatrix, and r^*i is the torque exerted by the inner disk on 
the librating region at the inner separatrix. We can write: 

and: 

= M\j{rp + Xs) - j{rp - Xs)] (6) 

where M is the mass flow from the outer disk to the inner disk, which would be 37rz/Soo in 
a strictly stationary regime, and where we have assumed: \dv/dr\ ^ \v/r\. To lowest order 
in Xs/rp, Eq. (6) can be recast as: 

Fr = MlQpTpXs + 2{rp + Xs)v+ - 2{rp - Xs)v_] (7) 

where v± is the perturbed azimuthal velocity at the outer (rcsp. inner) separatrix. 
Therefore, the use of Eqs. (4), (5) and (7) enables one to estimate the corotation torque 

has described at least A'" — 1 orbits, and at most + 1 orbits; taking N —>■ oo, one is led 
to (Q) = flp, where (Q) is the average angular velocity of the fluid element). In a viscous 
case, the positive torque arising from the close encounters of the outer disk fluid elements 
with the planet as they get sent to the inner disk needs to be taken into account, even if 
these fluid elements are only temporarily corotating with the planet. Failing to do so would 
in particular make it impossible to reconcile the torque expression given in this paper with 
previous estimates, since this positive torque almost cancels out the negative viscous torque 
exerted on the librating region; the net torque, which we identify as the coorbital corotation 
torque, comes from a balance between these two. 



3^ fdv\ 
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only from the knowledge of the unperturbed disk characteristics and the flow properties at 
the separatrix, provided a steady flow solution is known, since the steady flow assumption 
is necessary to write Eq. (3). 



5. Steady state solution 

In order to find a steady state solution, we write a radial viscous iffusion equation for 
the surface density. This equation needs to include the turnover of the horseshoe orbits. 
One simple way to do that is to include source and sink terms for \x\ < Xg- This is detailed 
in Appendix A. 

When this steady flow solution is used to evaluate the torque expression given by 
Eqs. (4), (5) et (7), in which we neglect the perturbed azimuthal velocity profile and the 
even part in x of the density profile perturbation (the 'dip' which surrounds the orbit), one 
obtains: 

Tc = \xlnp:^J^{zs) (8) 

where: 

where is a linear combination of the Airy functions Ai and Bi defined by Eq. (A18) and 
where: 

In particular, in the highly viscous case {zs -C 1), we have: 

^ „4r>2 



rc = -xtn;j:^ (u) 
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where we use the fact that (Abramowitz & Stegun, 1972): 

Eq. (11), with our notations, is the same as the expression given by Ward (1992): 

r'c = Ix^j:.. - (13) 

in which we have to set the slope /? = — dlogE/dlogr to since we are in the very viscous 
regime in a uniform surface density disk. 

In the general case, at lower viscosities, there is an ambiguity in the choice of the slope 
of P that has to be used in Eq. (13). Indeed, if we choose: 

A = atx = (14) 



dx 



then we get the expression: 



whereas if we choose: 



(15) 



^^^_E(x,)-E(-ag (16) 



2,Xg ^oo 



then we get: 



The graphs of Tc, and T'q^ are represented in Fig. 2. Tq^ is always smaller than Tc, 
and F'q^ is always larger than Fc- This means that it is always possible to find a value for 
(3 between /3i and /?2 for which Eqs. (8) and (13) give the same value. The ambiguity in 
the choice of /3 for Eq. (13) comes on one hand from the fact that this equation has been 
derived by an integration of the torque on all the librating streamlines of the horseshoe 
region, and on the other hand from the fact that the solution of Eq. (A17) does not lead 
to a constant slope over the whole horseshoe region. Such an ambiguity does not exist in 
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10.0000 



Fig. 2. — Coorbital corotation torque as a function of viscosity. The solid curve is a plot of 
expression (8), while the dashed and dotted ones represent respectively Eqs. (15) and (17). 
Here T^"^ stands for ^nlxjE^. 
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Eq. (8) since this latter has been derived taking into account only the disk perturbation at 
the separatrix. What happens in the interior of the librating region does not need to be 
known, provided we are dealing with a steady flow. 

The ratio TZ — i'rp/{VLpx1) = l/{27rZg) is also ^7^, where: 

^ Hlf ^^^^ 



and: 



'Tvisc (■^^) 



are respectively the turnover time of the outermost horseshoe orbit and the viscous diffusion 
time across the horseshoe region. From Fig. 2 we see that the coorbital corotation torque 
reaches half of its maximum value at TZ — v/x^g ~ 0.05, that is to say for ths ~ ^visc- When 
the turnover time is much larger than the viscous time, i.e. when TZ ^ 0.05, the corotation 
torque reaches its maximal value, since libration is inefficient at removing the vortensity 
gradient. On the contrary, when the viscous time is much larger than the turnover time 
(i.e., when TZ <^ 0.05) libration can remove the vortensity gradient and the corotation 
torque has a negligible value. 



5.1. Cut-ofF at high viscosity 

As discussed in section 3.2, the outer separatrix can be approximated as a branch of 
parabola in the 6*, r plane, and so far we have assumed that the viscosity is small enough 
so that we can approximate the separatrix by a line of constant distance to the orbit. 
However, if the viscosity is very high this is no longer true, and the separatrix can reach the 
orbit before having described a complete revolution in the planet frame. From Eq. (2), we 
can give an estimate of the critical value for the viscosity above which this will happen, by 
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setting 9o — 0, X — Xs and 6 — 27r. This leads to: 

= ^ (20) 

The 2D analysis we have developed here is mainly valid for mildly embedded protoplanets, 
for which the Hill radius is a sizable amount of the disk thickness; this would correspond 
to values of the a-parameter of about a ~ 0.1, probably too high for protoplanetary disks, 
and therefore this regime is not very relevant for the analysis described in this paper. 



The coorbital corotation torque general expression in a steady regime, which we have 
given in Eqs. (4), (5) and (7), contains additional terms that we have so far neglected. In 
the following we still neglect the terms linked to the perturbed azimuthal velocity v, the 
analysis of which will be found in Appendix B, but we do not neglect, as we have done 
before, the perturbation of the density profile in the vicinity of the planet orbit. Indeed, 
as is seen frequently in numerical simulations, the mildly embedded planets, which are not 
massive enough to fulfill the gap opening criteria, are nevertheless surrounded by a shallow 
axisymmetric depletion of their coorbital region. As a consequence, a new term arises in 
the torque expression, which reads: 



where Eg is the average of the surface density at the outer and inner separatrices (in 
different words it is the value of the even part in x of the perturbed density, sampled at any 

scparatrix). If we assume that this value is the same as the value at the orbit (which means 
that the dip edges lie further from the orbit than the separatrix), then we can write: 



6. Additional terms and Lindblad torque coupling 



= 37ri/(Eoo — Es 




(21) 




(22) 
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where r^ij > is the one-sided Lindblad torque. Eq. (22) is obtained by equating the 
viscous torque and the torque given by the wake damping, and integrating it between the 
orbit and an outer radius where most of the wave angular momentum has been transfered 
to the disk (i.e., it is obtained by integrating twice the even part in x of Eq. (A9) in the 
steady state case). Because of this integration, the assumption of local damping is not 
necessary in writing Eq. (22). However, the damping must occur over a radial range smaller 
than the planet orbital radius for this expression to be accurate. Recent work by Goodman 
& Rafikov (2000) shows that this may be indeed the case even for relatively small planet 
masses. Using Eq. (22), the torque expression of Eq. (21) can be recast as: 

= -Vlr (23) 

It can be argued that this result has been obtained assuming a strictly steady flow, and in 
particular that the mass flow from the outer disk to the inner disk across the planet orbit 
is exactly Sttz/Eoq. A higher mass rate than this value (which is 37ri/Eoo either far from the 
planet or in the unperturbed disk) would lead to a value for Vq higher than the one given 
by Eq. (23). The analytical work for an inviscid disk by Lubow (1990), and Goodman & 
Ravikov (2000), as well as numerical evaluation of the mass flow rate accross the orbit of 
a non-accreting protoplanet (Masset and Snellgrove, 2001, in their Fig. 4) suggest that 
the actual mass flow rate across corotation is higher than 37ri/Eoo, and that the difference 
between the actual value and 37rz/Soo scales with HooMp. This issue definitely requires 
further work, as it could provide another positive additional term to the corotation torque. 

The other terms which appear in the general torque expression are evaluated in B. 
They all correspond to the difference between the outer and the inner separatrix of 
quantities involving the perturbed azimuthal velocity v. The perturbed azimuthal velocity 
can be split in two parts: the part which corresponds to the disk axisymmetric profile 
perturbation under the action of the Lindblad torques, which is odd in a; to lowest order in 
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the disk aspect ratio, and an even part in x (to lowest order in Xs), which corresponds to the 
odd perturbation density due to hbration that we have analyzed at Appendix A. However 
in the analysis of Appendix A we have neglected the perturbation of the azimuthal velocity, 
and therefore we cannot estimate the additional terms linked to the even component of the 
perturbed azimuthal velocity. Taking this latter into account would have led to a higher 
order diffusion equation in Eq. (All) which would have involved the disk thickness, which 
has not been considered in the simple analysis of Appendix 5. This will be done in a 
forthcoming work. We note that the additional term linked to this even component of the 
perturbed azimuthal velocity may play a non negligible role and significantly modify the 
viscosity dependent part of the corotation torque given by Eq. (8). 

The additional term linked to the odd component of the perturbed azimuthal velocity 
depends only on the value of this latter at the separatrix. However, one can only give a 
rough estimate of the maximal value of these terms. Their actual value will depend on the 
exact position at which the separatrix samples the perturbed velocity profile. It turns out 
that the maximal value of this term is of the same order of magnitude as the additional 
term of Eq. (23): 

V'S<\-{'i-mLR (24) 

Hence we see that the main effect of the additional terms (both Tq and Tq) is to introduce 
a coupling with the one-sided Lindblad torque, as they add a small, positive fraction of this 
latter to the coorbital corotation torque. Now, the differential Lindblad torque is a small, 
negative fraction of the one-sided Lindblad torque, which scales with the aspect ratio (for 
instance in the disk we consider, with a uniform surface density and a constant aspect 
ratio h' , ^Tlr ^ —ShTLR)- Therefore this additional term plays against the differential 
Lindblad torque and one has to expect that the migration is considerably slower than 
standard type I migration estimates, or even reversed, over a certain planet mass range, 
in very thin viscous disks. The exact determination of this mass range has to be achieved 
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through numerical simulations, since the analytical determination of the exact perturbed 
density profile, of the exact azimuthal perturbed velocity profile, and of the exact separatrix 
position seems a hardly tractable task. 



We have derived the expression of the coorbital corotation torque exerted by a uniform 
viscosity and surface density disk on an embedded protoplanet on a fixed circular orbit, 
assuming a 2D steady flow in the planet frame. We have obtained the following expression 
for the torque: 



where J-{zs) is given by Eq. (9), and where Q{xs) is a function of the order of unity, 
the exact value of which depends on the position of the separatrix and of the perturbed 
axisymmetric profile of the disk under the action of the Lindblad torques, and the maximum 
value of which is 7/3, as can be inferred from Eqs. (23) and (24). The first term in the 
R.H.S. of Eq. (25) corresponds to previous estimates of the corotation torque, and exhibits 
an explicit dependence on the viscosity: when the viscous diffusion timescale across the 
coorbital region is larger than the horseshoe libration time, this term is negligible, while 
when the viscous diffusion time is shorter than the libration time, it reaches its maximum 
value. The second term of the R.H.S. of Eq. (25) comes from a dependence of the viscous 
torque acting on the boundary (separatrix) of the trapped librating region on the perturbed 
axisymmetric profile of the disk under the action of the one-sided Lindblad torques on 
each side of the orbit. This term slows down the migration and could even reverse it in 
sufficiently thin disks. 



7. 



Discussion 




(25) 
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This analysis is valid as long as the flow around the planet can be considered as a 
2D steady flow. Therefore the planet Hill radius needs to be a sizable amount of the disk 
thickness, and gas accretion cannot occur onto the planet (otherwise the flow cannot be 
considered as steady if the accretion time is short, and the separatrix splits into two lines 
which enclose the band of the material which is going to be accreted; see e.g., Lubow et 
al. 1999). Therefore the mass limit for this analysis is roughly 15 M^. Another condition 
for the analysis presented above to be valid is that the planet does not open a gap in the 
disk (otherwise the corotation torque is negligible). The depletion which appears around 
the planet orbit and which leads to the Lindblad torque coupling needs to be shallow. 

The analysis presented here rehes on the existence of a viscous drift of the disk 
material across the planet orbit. Therefore it can apply to those parts of the disk where 
MHD turbulence occurs and leads to an effective viscosity (Balbus & Hawley, 1991). 
Hence it is likely to apply mainly to the innermost part of the protoplanetary disks (the 
central astronomical unit). It is questionable whether this kind of turbulent transport can 
be correctly modeled by an effective viscosity when applied to the disk-planet angular 
momentum exchange in the coorbital region. In particular the steady state assumption that 
was used to estimate the corotation torque needs to be reconsidered. 

One may speculate about the extrapolation of the calculations exposed here to smaller 
masses, for which the Hill radius is much smaller than the disk thickness (deeply embedded 
planets). However the 2D flow assumption was essential for the identification of a closed 
region of trapped material. There is no guarantee that we could perform the same kind of 
calculation in the 3D case. However, if the planet has no inclination and no eccentricity, 
the gas motion in the disk equatorial plane would still be 2D for symmetry reasons, and 
therefore would have a set of closed streamlines in the steady state regime. However, the 
extension of the trapped region above the equatorial plane would need to be discussed. 
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In the case of a finite excentricity, the analysis exposed here is no longer valid; in 
particular, the planet and the disk can exchange angular momentum at the excentric 
corotation resonance locations, which are not localized where the disk material corotates 
with the planet guiding center (i.e., which do not lie at the orbit semi-major axis, if one 
neglects the radial pressure gradient as a partial support to gravity in the disk). However, it 
would be interesting to find a link between the torque expression obtained by summing over 
the resonances and an analysis similar to the one exposed above, in which the separatrix 
would become a band of width 2ea (using the fact that the epicyclic frequency is much 
higher than the outermost horseshoe turnover frequency). In particular it would be of 
interest to see how the coupling with the one-sided Lindblad torque is modified by the 
introduction of a finite excentricity. 

This analysis needs to be extended to the case of an arbitrary surface density profile, 
and to a migrating planet (for which the mass flow across the orbit will therefore differ from 
37ri/Eoo). In this case the corotation torque would include a delayed term proportional to d, 
the planet drift rate, and would lead to a feedback on the migration. This will be presented 
in a forthcoming work. 

A. Steady flow solution 

In the horseshoe region we write the continuity equation integrated over an arc of a 
ring [r, r -|- Sr] which excludes the horseshoe U-turn ends of the streamhnes. We are led to: 

f-^^-^(7l"-"<-|..-.l.)=o 

where u is the radial velocity. In Eq. (Al) we assume that the arc of ring spans almost an 
angle 2%, and that it is mapped, after a horseshoe U-turn, onto the arc of ring of radius r' 
and width Sr'. We also neglect the perturbed azimuthal velocity. The variables S' and fl' 
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stand respectively for S(r') and Q{r'). The left term of the bracket corresponds to the mass 
which comes from 'the other side' (i.e., the ring [r',r' + Sr']), and enters the ring [r, r + Sr] 
after a horseshoe U-turn, while the second term in this bracket corresponds to the mass 
which leaves the ring [r, r + 5r]. 

The torque per unit mass T acting on the disk material can be decomposed as the sum 
of the viscous torque per unit mass and the torque per unit mass {27r RE)~^ dT / dr arising 
from the wake damping: 

1 d{r^uE§) 1 dr 

rE dr 2nREdr ^ ' 

Writing T — udj/dr — udjo/d — urQ/2, and using Eqs. (Al) and (A2), one is led to the 
following radial diffusion equation: 

^ = ^('-\n -n'lE'^-ifi -i^is'Vf^Vs-^-^— (A3) 

dt 27r \ r ^ 5r ^ J \Xs J dr'^ nr'^fl dr"^ 

where 11 is defined as: 

n : u I — > 1 if \u\ < 1 

otherwise 

and accounts for the fact that the source and sink terms are present only over the 
horseshoe region (i.e. \x\ < Xg). In Eq. (A3) we have assumed r|9^E/9r^| » \dT,/dr\ and 
r\d'^T/dr'^\ \dT/dr\. If we neglect both pressure and viscous effects during a U-turn at 
one end of the horseshoe orbit, then the conservation of the Jacobi constant of a fluid 
element implies, as has been shown by Ward (1992), that: 

Sr r'B'\n'-nJ 



5r' rB\n-np\ 



(A4) 



where B, the second Oort's constant, in the case of a Keplerian disk, is also B — Vt/A. 
Using Eq. (A3) and ( A4) , one can write: 

-T- = — S' E n„ - Q n — + 3z/— (A5) 

dt 27r V ^' ^ \xs dr^ nrpp dr^ ^ ' 
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If we write x' — r' — Vp, and if we define the dimensionless distances to tfie orbit as x — x/vp 
and x' — x'/fp, then we have: 

§ = 1 - + ^x' + + \x- + 0{x') + 0{x'') (A6) 

From Eq. (A4), we can infer that x — —x' to lowest order in £, therefore: 

where C is a numerical constant. Hence the bracket in the source term of Eq. (A5) reads: 
E'^-E = E(-x) - E(x) - 3xE(-x) (A8) 

1 Fit 

+ -(9 + 6C)£'E(-i;) + Cx^-qt{-x) + 0(i;3) 
If we keep only terms to first order in x in this expression, Eq. (A5) can be recast as: 

— = -n,[E(-x)-E(x)-3xE(-x)] (A9) 

We call E and E respectively the odd and even part in x of E: 

t{x) = i[E(x) -E(-x)] 
t{x) = |[E(x) + E(-x)] 

Wc write E(a;) = Ec + 0(a;^), where Ec = E(0) = E(0). Taking the odd part in x of both 
sides of Eq. (A9), and keeping terms only up to O(x^), one is led to: 

where we have neglected, at this step, the Inner and Outer Lindblad torque imbalance, and 
thus the odd part in a; of F and its second derivative vanishes. We write s = E/Ec and 
i> = ^/{rpflp), and we assume that we have a steady state. Therefore, for x > 0, we have: 

d'^s 1 

z>— - = — [2s + 3£1£ ifx<Xs (All) 
ox^ 47r 

= otherwise 
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The first line of Eq. (All) shows the balance between hbration, which tends to create a 
profile of density such that s = —fa;, and the viscous diffusion. We want to find a solution 
to this ODE valid for any £ > 0, and such that s(0) = 0, since s is an odd function of x. 

The general solution of Eq. (All) for x > Xg is s = ax + bx'^, which does not diverge 
for large values of x only if a = and 6 = 0. Therefore, we are led to seek a solution of 
Eq. (All) for x e [0,Xs] which fulfills the two conditions: 

s(0) = 
s'{xs) = 



where s' = ds/dx. If we define 

■9 = .s + 

2 



where: 



(A12) 



S = s+-x (A13) 



and 

z = (27ri>)-^/2£, (A14) 
then Eq. (All) can be recast, ior < z < Zg — {27rC')~^/^Xs, as: 

^g^^zS (A15) 

and the set of conditions (A12) becomes: 

S{0) = 
§{zs) = |(27r^)l 

The solution of Eq. (A15) which fulfills the conditions (A16) is: 



(A16) 



, , Bi(^) - VSAi(z) 
g(z) = ^ , ^ ^ (A18 

^ 2Bi'(0) ^ ^ 

is the linear combination of Airy functions Ai and Bi such that g{0) = and g'{0) = 1. This 
solution can be used to evaluate the coorbital corotation torque given by the expressions (4) , 
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(5) and (7). In this part we neglect the perturbed azimuthal velocity, and we neglect as well 
the even part of the density profile perturbation (the shallow 'dip') which appears around 
the orbit for a sufficient mass (i.e., we write E = T^oo)- Furthermore we assume that the 
mass flow across the coorbital region from the outer disk to the inner disk is M = Stti/'Eoo- 
Therefore the expression of the corotation torque, in this approximation, reduces to: 

Tc = 37ru[(rp-Xsfn{rp + Xs)^J:s+ (A19) 
- (rp - Xsf^{rp - Xs)T.S'] 

which leads to: 

Tc = -67rz/rJOpEoos(£,) (A20) 
Using Eqs. (A13), (A17) and (A20), one can flnally deduce Eq. (8). 



B. Additional terms 

The additional term arising from the even part in x of dv/dr in Eq. (5) reads: 

r'c' = AnuTpX^J:, i^—j (2 - (3) (Bl) 

where Eg is the surface density at the planet orbit. 

In order to give an estimate oi v, we write the axisymmetric part of the linearized 
radial Euler equation, in which we omit the gradient of the m — perturbed potential of 
the planet. This reads: 

2nv-cl^^0 (B2) 

where is the sound speed. We can give a maximal estimate of ^^E by writing that the 
wave damping occurs over a radial range of at least one disk thickness. This reads: 

3-S^'-^f < (B3) 
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Using Eqs. (B2) and (B3) one infers: 



V < „ ^ Vlr (B4) 



In the same manner one can write an upper limit for the perturbed azimuthal velocity 
derivative: d^v <v/h. Therefore: 

dr ~ 67ri/Eoor^ ^ ^ 

Assuming that Eg ~ E^o, i.e., the dip surrounding the planet orbit is shallow, one can infer 
from Eqs. (Bl) and (B5) that: 

r^'<^(2-/3)-r^H (B6) 

The additional term arising from v in Eq. (6) reads: 

r^^^ = 37ri/Eoo • 2x^1; (B7) 
Using Eq. (B4), one finds as a maximal value for this term: 



r'c < ^Vlr (B8) 



Since h <^ rp and Xs <^ Tp, this term is negligible. 
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